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We present a study of the P-wave spin -triplet charmonium XcJ decays (J = 0, 1, 2) into 7r°7r 7r 7r°. 
The analysis is based on 106 million ip' decays recorded with the BESIII detector at the BEPCII 
electron positron collider. The decay into the 7r°7r 7r 7r° hadronic final state is observed for the first 
time. We measure the branching fractions B(xco ttWV ) = (3.34±0.06±0.44) x 10" 3 , B(xa 
ttWV ) = (0.57 ± 0.03 ± 0.08) x 10" 3 , and P(xc2 -> ttVW) = (1.21 ± 0.05 ± 0.16) x 10" 3 , 
where the uncertainties are statistical and systematical, respectively. 

PACS numbers: 13.25.Gv, 14.40.Pq, 13.20.Gd 



I. INTRODUCTION 

In the quark model, the \cJ (J = 0, 1, 2) mesons are 
the 3 Pj charmonium states. Their decays are experimen- 
tally and theoretically not as well studied as the vector 
charmonium states J/ip and ip' . In contrast to the latter 
ones, XcJ cannot be produced directly in e + e~ annihila- 
tion. However, they can be produced in radiative decays 
■ip' — > "/XcJ j providing a clean environment to study their 
decays. 

Recent theoretical work indicates that the Color Octet 
Mechanism [l[ could have large contributions to the de- 
cays of the P-wave charmonium states. However, these 
calculations as well as experimental measurements still 
have large errors and thus more precise experimental data 
besides more theoretical efforts are mandatory to further 
understand XcJ decay dynamics. Furthermore, the XcO 
and Xc2 states are expected to annihilate via two-gluon 
processes into light hadrons and may therefore allow the 
study of glueball dynamics. Thus the measurement of as 
many exclusive hadronic XcJ decays as possible is valu- 
able. 

The XcJ decays into four pions have the largest branch- 
ing fractions among the known hadronic XcJ decay modes 
Presently only the decays into it + tt~7t + it~ and into 
7r + 7r~7r°7r are measured by previous experiments. The 
branching fractions are shown in Table [I] In this paper, 
we present a study of exclusive XcJ decays into 7r 7r 7r 7r . 



channel 




branching fraction [%] 


XcO — > TT + TT~ 


~TT + TT~ 


2.27 ±0.19 


Xcl — > TT + TT~ 


~TT + TT~ 


0.76 ±0.26 


Xc2 -> TT + TT~ 


~7T + 7r~ 


1.11 ± 0.11 


XcO -> TT + TT~ 


-JJ 
7T 7T 


3.4 ±0.4 


Xcl — > TT + TT~ 


TT TT 


1.26 ±0.17 


Xc2 -> TT + TT~ 


TT TT 


2.00 ± 0.26 



TABLE I: Branching fractions of x»J into tt + tt tt + tt and 



II. THE BESIII EXPERIMENT AND DATA SET 

We use a data sample of about 106 million ip' decays 
recorded with the BESIII detector 3] at the energy- 
symmetric double-ring e + e~ collider BEPCII [4j. The 
primary data sample corresponds to an integrated lumi- 
nosity of 156.4pb _1 collected at the peak of the ip' res- 
onance. In addition, a 42.6 pb" 1 data sample collected 
about 36 MeV below the resonance is used for background 
studies. 

The BESIII detector is described in detail elsewhere 
Charged particle momenta are measured with a 
small-celled, helium gas-based main drift chamber with 
43 layers operating within the 1 T magnetic field of a 
solenoidal superconducting magnet. Charged particle 
identification is provided by measurements of the specific 
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ionization energy loss dE/dx in the tracking device and 
by means of a plastic scintillator Time of Flight system 
composed of a barrel part and two end caps. Photons are 
detected and their energies and positions measured with 
an electromagnetic calorimeter (EMC) consisting of 6240 
CsI(Tl) crystals arranged in a barrel and two end caps. 
The return yoke of the magnet is instrumented with Re- 
sistive Plate Chambers arranged in 9 (barrel) and 8 lay- 
ers (end caps) for discrimination of muons and charged 
hadrons. 



III. DATA SELECTION 

We reconstruct the entire event from the decay chain 
of the charmonium transitions ip' — > 7Xcj followed by the 
hadronic decays Xcj —> 7r 7r°7r 7r°. A photon candidate is 
defined as a shower detected with the EMC exceeding an 
energy deposit of 25MeV in the barrel region (covering 
the region | cos 9\ < 0.8 of the polar angle) and of 50 MeV 
in the end caps (0.86 < |cos6>| < 0.92). The average 
event vertex of each run is assumed as the origin of these 
candidates. We restrict the analysis to events having nine 
photon candidates and no reconstructed charged particle. 
The energy sum of the nine photons must be within the 
range 3.45-3.80 GeV. To reconstruct 7r° — > 77 candidates 
we use pairs of photon candidates, having an invariant 
mass between 110 and 150MeV/c 2 . Fig. [T] shows the in- 
variant 77 mass distribution of all photon pair combina- 
tions for selected 97 events. A clear ir° signal is visible. 
By combining four ir° candidates with an additional pho- 
ton candidate being detected in the EMC barrel, where 
the same photon candidate must not be used more than 
once, the complete event is reconstructed. Different pair- 
ings of the photon candidates can yield more than one 
77T 7r 7r 7r° candidate per event. Therefore we use the 
pairing which leads to the minimal 

2 _ \ " Kyj ~ TO^o) 2 

i 

calculated from the invariant mass m 77i i of the i th tt° 
candidate for a given 77r 7r°7r 7r° candidate, the nomi- 
nal 7T° mass Q m n o, and the ir° — > 77 invariant mass 
resolution a of 6.5 MeV/c 2 . Combinatorial background is 
suppressed strongly by demanding jq^ < 15. 

Potential backgrounds can arise from the transition 
tp' — > -K 0/ K J/ij) followed by hadronic or radiative decays 
of the J/0 to final states with higher photon multiplicity. 
We therefore require the recoil mass wir of any di-7r° pair 
with respect to the ip' to be \mn — mju, \ > 100 MeV/ c 2 , 
where mj/^ is the nominal J/ip mass |5(. 

The spectrum of the energy E* of the photon from 
the ip' radiative transition in the center of mass frame 
is shown in Fig. [5] Clear x c0 , x c i,and Xc2 signals with 
small background are evident. Analysis of the continuum 
data sample yields only four events passing the selection 




m 7Y [GeV/c 2 ] 



FIG. 1: Invariant 77 mass distribution of all photon pair 
combinations for selected 97 events. The arrows indicate the 
mass window used for the selection of n° candidates. 




E* y [GeV] 

FIG. 2: The spectrum of the energy E* f of the radiative pho- 
ton from ip' — y 7Xcj with the result of the fit (solid curve) 
described in the text. The dashed curve shows the back- 
ground line shape and the dotted curves represent the signal 
line shapes derived from MC simulations (see Sect. IV)) . 



criteria and thus reveals no significant background. Peak- 
ing components of the background are investigated from 
simulated Monte Carlo (MC) events and are discussed in 
Sec. ED 

We further look into resonant substructures in the 
7r 7r°7r 7r° final state. The production of intermediate 
resonances in the decay could have an impact on the de- 
tection efficiencies. Here, we only investigate the gross 
substructures by plotting the invariant mass of any di-7r° 
pair in the final state versus the corresponding mass of 
the other pair for the three Xcj signal regions (Fig. [3]). 
The defined x c0 , % cl , and Xc2 signal regions correspond 
to E* energy ranges of 220-290 MeV, 160-180 MeV, and 
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115-135 MeV, respectively. As seen in Fig. [3] production 
of Kg and /o(980) in the Xco decay is evident. Accu- 
mulation of events is also observed in the mass region 
around 1300, 1700, and 1950 MeV/c 2 . Structures around 
1300MeV/c 2 are observed in Xci and Xc2 decays. All 
these structures need further careful investigation using 
partial wave analysis techniques with increased data sam- 
ples being collected in the future. As for Xco decays, 
production is also observed in Xc2 decays, while the decay 
Xci —> KgKg is forbidden by parity-conservation. 

For the measurement of branching fractions we in- 
clude all sub-resonant decay modes but explicitly exclude 
the Xco and Xc2 decay mode to K^K S . Therefore we 
reject events where the invariant mass m\i of any di- 
7T° pair and the invariant mass 77134 of the correspond- 
ing other di-7r° pair of the 77r°7r 7r 7r° final state fulfills 

\ {m\2 — m K o) 2 + (77134 — m K o) 2 < 100 MeV/ c 2 , where 

L is the nominal K s mass fi. 



IV. MONTE CARLO STUDIES 

A detailed MC simulation of the BESIII detector based 
on GEANT4 [6] is used to determine efficiencies, signal 
shapes, and background contributions. The production 
of the ifj' resonance is simulated using the kkmc event 
generator Q. Decays of the 7// and subsequent parti- 
cles in the event are modeled by evtgen Q. Simulated 
events pass the same reconstruction algorithms and se- 
lection criteria as data. 

Signal MC data samples of 500k events for each de- 
cay tp' —> JXcJ, XcJ 7T 7T 7T°7T are generated using a 
1+A cos 2 (#) distribution, where 8 is the angle between the 
direction of the radiative photon and the positron beam, 
and A = 1, —1/3, 1/13 for J = 0, 1,2 in accord with ex- 
pectations for electric dipole (El) transitions. The XcJ 
decay products are generated using a flat angular distri- 
bution. Intrinsic width and mass values as given in Q are 
used for the XcJ states in the simulation. The obtained 
efficiencies for Xco, Xci, and Xc2 are (10.16 ± 0.05)%, 
(11.54 ± 0.05)%, and (10.85 ± 0.05)%, respectively, in- 
cluding detector acceptance as well as reconstruction and 
selection efficiencies. 

In addition we use MC data samples to investigate 
sources of the peaking backgrounds. For each of the stud- 
ied XcJ decay modes listed in Table [TT] we generated at 
least 100k events. The contribution of the total peaking 
background is estimated from a fit to the reconstructed 
E* spectrum. The fit procedure is the same as applied 
for data and will be addressed in Sect. [Vj The largest 
peaking background contribution is found to come from 
Xci 777r°7T° decays, with 77 — > 7r 7r°7T°, where one of the 
7T° has low momentum and is not detected. 



channel 




n XaO 




n Xc2 


XcO - 


>K%K%K%- 


-> TT TV 


1.6 


0.3 





XcO - 


► VI 




0.2 








Xci - 






1.2 


45.2 





Xci - 


> -yJ/tp , J/0 - 








1.3 





Xci - 


> fJ/ip , J/0 - 










0.1 


Xci - 


>K S K S ,K S - 










0.6 


Xci - 


► 777rV° 










3.8 


Xci - 


> 7777 













Xci - 


> 7 J/0 , J/0 - 










0.6 


Xci - 


> 7 J/0 , J/0 - 













Sum 






3.0 


46.8 


5.1 



TABLE II: Expected number of background events peaking 
at the XcJ signal regions as derived from MC simulations. 



V. FITTING PROCEDURE AND EXTRACTION 
OF BRANCHING FRACTIONS 

The E* spectrum shown in Fig. [2] is fitted using an 
unbinned maximum likelihood fit. The XcJ signal line 
shapes are extracted from the MC simulation. A 2 nd 
order Chebychev polynomial is used to describe the non- 
peaking background. From the fit Xco, Xci, and Xc2 signal 
yields of 3299 ± 67, 655 ± 32, and 1169 ± 41, respectively, 
are obtained. 

To access the goodness of the fit, we repeat the fit to 
the E* spectrum using a binned least-squared fit. Ap- 
plying a binning of 5 MeV, this fit yields a x 2 value of 
102 with 84 degrees of freedom. 

The fit does not account for the peaking component 
of the background. We estimate the number of events 
peaking at the position of the XcJ signals by fitting the 
E* spectrum derived from the MC samples generated 
for background studies. The same fitting procedure as 
for data is applied, except that the parameters of the 
polynomial function describing the shape of the non- 
peaking combinatorial background are fixed to the values 
obtained from the fit to data. From the extracted signal 
yields the expected number of peaking background events 
is calculated using the branching fractions of the XcJ de- 
cays as given in [5J . We relate the unmeasured branching 
fractions of the decays XcJ — ► t]it it to the corresponding 
branching fractions of XcJ — > i]tt + tt~ decays by the iso- 
spin ratio of the two final states. The estimated number 
of peaking background events for the investigated chan- 
nels is given in Table [TTT In total a peaking background of 
3.0, 46.8 and 5.1 events to the Xci), Xci, and Xc2 signals, 
respectively, is derived. The largest contribution arises 
from the decay Xci 7771° 7r°, amounting to 45 events at 
the Xci signal region. 

Although Xco -> K° s K a s and % c2 -4 K° S K° S decays for 
the branching fraction measurement were excluded, we 
considered these channels as a potential peaking back- 
ground source. The feed-through of the two decays to 
the XcO; Xci, and Xc2 signals is expected to be 1.6, 0.3, 
and 0.6 events, respectively. 
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FIG. 3: Shown is the invariant di-7r° mass plotted versus the other invariant di-7r° mass for events selected from the a) \c0, b) 
Xci, and c) \c2 signal regions. All possible ty° ty° combinations of the 7r 7i- 7r°7r hadronic final state are plotted. In addition 
the plots are symmetrized; thus each event enters six times to the plots. The one-dimensional projections are shown in d) e) 
and f) for \c0, Xci, and \ C 2, respectively, where the dots with error bars show the data and the shaded histograms show the 
di-7r° mass distributions obtained from signal MC events simulated without intermediate resonances for the \cj decays. 



The expected number of peaking background events is 
subtracted from the yields observed for data. These cor- 
rected yields N are then converted to branching fractions 
using 



} e • N$> ■ B{i>' -> 7Xc./) • -> 77 ) 4 (2) 



where e is the efficiency; NL is the number of ip' in the 
data sample; and B{\jj' — > 7 Xcj) and B(n a — > 77) are the 
branching fractions of radiative tj/ transitions into Xcj 
and of the decay ir° — > 77 [jj, respectively. The number 
of ip 1 and its combined statistical and systematical uncer- 
tainties are determined to be Nm = (1.06 ± 0.04) x 10 8 
1- 



VI. ESTIMATION OF SYSTEMATIC 
UNCERTAINTIES 

Several sources of systematic uncertainties are consid- 
ered for the measurement of the branching fractions, in- 
cluding uncertainties on the photon detection and recon- 
struction; the event selection; the fitting procedure and 
peaking background subtraction; and the number of tp' 
decays in the data sample. The investigated uncertain- 
ties are summarized in Table IIIII and will be discussed in 
detail in the following. 

a. Photon Detection The uncertainty due to pho- 
ton detection and conversion is 1% per photon. This is 
determined from studies of photon detection in well un- 
derstood decays such as J/ip — > p°ir° and the study of 
photon conversion in the process e + e~ — > 77. 

b. Event Selection By varying the requirement on 
xfvr, the 7T° mass window and the total energy of the 
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Xoo[%] 


Xd{%] 


Xta[%] 


photon detection 


9.0 


9.0 


9.0 


decay model 


6.3 


6.3 


6.3 


B(ip' -> 7Xcj) 
number of ip' events 


3.2 


4.3 


4.0 


4.0 


4.0 


4.0 


total energy 


3.0 


3.0 


3.0 




2.5 


2.5 


2.5 


reconstructed ir° mass 


1.4 


1.4 


1.4 


fitting range 


1.4 


2.9 


0.9 


signal line shape (energy resolution) 


0.5 


1.2 


0.4 


signal line shape (energy shift) 


0.1 


0.6 


1.0 


background shape 


1.0 


1.0 


1.0 


peaking background subtraction 


0.1 


0.8 


0.5 


MC statistics 


0.5 


0.5 


0.5 


trigger efficiency 


<0.1 


<0.1 


<0.1 


total uncertainty 


13.2 


13.6 


13.2 



TABLE III: Summary of the systematic uncertainties. 



77r°7r 7r 7r° candidates used for the event selection in 
data and MC events, we investigate the systematic un- 
certainties in modeling the distribution of these param- 
eters. The largest deviation of the branching fractions 
from the default values sets the scale of our systematic 
uncertainty, and we assign a uncertainty of 2.5% for the 
xl-xi 3% for the total energy; and 1.4% for the ir° mass 
window requirement. 

c. Monte Carlo Decay Model The efficiencies for the 
processes ip' — > "iXcj, Xcj — > 7r 7r°7r 7r° are determined 
from MC simulations, where no intermediate resonances 
and a flat angular distribution have been considered in 
the Xcj decays. As discussed in Sect. |TTT] this analy- 
sis reveals the presence of intermediate resonances in 
the Xcj 7r°7r 7r 7r° decays. This could have an im- 
pact on the detection efficiencies, which we consider as 
a systematic uncertainty. Wc determine the efficiencies 
from our simulations including the sub-resonant modes 
/o(980)/ (980), / 2 (1270)/ 2 (1270), / (1370)/ (1500), and 
/o(1370)/ (1710) for X c0- We considered / 2 (1270)7r°7r° 
for Xd and / 2 (1270)/ 2 (1270) and / (1370)/ (1710) for 
Xc2- We find there is no large efficiency difference 
from that of phase space. The largest difference with 
respect to the efficiency obtained for the simulation 
without intermediate resonances is observed for Xco ~ > 
/o(1370)/o(1710) decay. We take this difference as a con- 
servative estimate of the uncertainty due to the MC decay 
model and assign a uncertainty of 6.3% for the Xco, Xci, 
and Xc2 branching fractions. 

d. Fitting Procedure The Xcj yields determined 
from the fit to the E* spectrum determines the branching 
fractions. We repeat the fit with appropriate modifica- 
tions to estimate the systematic uncertainties due to the 
fitting procedure. The difference of the derived branch- 
ing fractions with respect to the values derived from the 
standard fit is considered as a systematic uncertainty. 

We smear the resolution function of the Xcj signals ob- 
tained from MC simulations by 1% to estimate the sys- 



tematic uncertainties of modeling the photon resolution 
in the MC simulation and shift the signal mean values by 
±1 MeV to estimate the systematic uncertainties due to 
the absolute energy calibration. The assigned uncertain- 
ties are given in Table Mil To estimate the uncertainty 
due to the non-peaking background parametrization we 
use a third order instead of a second order Chebychev 
polynomial. This uncertainty is found to be 1%. For the 
nominal fit, the E* spectrum is fitted in the interval 0.05- 
0.5 GeV. A series of fits using different E* intervals is per- 
formed and the largest change of the individual branch- 
ing fractions is assigned as a systematic uncertainty (see 
Table HD. 

e. Peaking Background Subtraction The number of 
peaking background events is estimated from MC simu- 
lations and is subtracted from the signal yields obtained 
from the nominal fit. The major source of peaking back- 
ground is the decay Xci — > i]ir°ir . To determine the 
number of background events, the branching fraction of 
this decay mode is computed from the Xd t]it + Tr~ 
branching fraction exploiting the isospin relation of the 
two decays. The uncertainty of the branching fraction of 
Xci r]TT TT Q leads to a systematic uncertainty of 0.8% 
for Xci ■ For Xco and Xc2 we take the number of subtracted 
background events as systematic uncertainties. 

/. Other Systematic Uncertainties For the normal- 
ization of the branching fractions, the number of ip' 
events in the data sample determined according to the 
method as given in is used. This method yields a sys- 
tematic uncertainty of 4%. The uncertainty due to the 
branching fractions ip' — > IXcj is 3.2% for Xco, 4.3% for 
Xci, and 4.0% for Xc2- 

A systematic uncertainty of 0.5% is assigned due to 
the statistical error of the efficiencies as determined from 
MC simulations. 

The systematic uncertainty due to the simulation of 
the trigger efficiency is found to be less than 0.1% (loj . 

g. Total Systematic Uncertainty We assume that all 
systematic uncertainties given above are independent and 
add them in quadrature to obtain the total systematic 
uncertainty. 



VII. CONCLUSION 

In summary we have measured the branching frac- 
tions of Xd -> 7r 7r°7r 7r° decays B(xc0 -> 7r°7r 7r 7r°) = 
(3.34 ± 0.06 ± 0.44) x lO" 3 , B(x c i -> 7r 7r°7r 7r°) = 
(0.57 ± 0.03 ± 0.08) x 10~ 3 , B( X c2 -4 7r°7r 7r 7r°) = 
(1.21 ± 0.05 ± 0.16) x 10~ 3 for the first time, where the 
quoted uncertainties are statistical and systematical, re- 
spectively. The 7r°7r°7r 7r hadronic final state contains 
a rich substructure of intermediate resonances. The re- 
ported branching fractions include decay modes with in- 
termediate resonances except Xco — > K^K® and Xc2 — > 
K®K®, which have been removed from this measurement. 
Our observation improves the existing knowledge of the 
Xcj states and provides further insight into their decay 
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mechanisms. Based on our results detailed studies of the 
sub-resonant decay structure with increased data samples 
may follow in the future. 

VIII. ACKNOWLEDGMENTS 

The BESIII collaboration thanks the staff of BEPCII 
and the computing center for their hard efforts. This 
work is supported in part by the Ministry of Science and 
Technology of China under Contract No. 2009CB825200; 
National Natural Science Foundation of China (NSFC) 
under Contracts Nos. 10625524, 10821063, 10825524, 
10835001, 10935007; the Chinese Academy of Sciences 



(CAS) Large-Scale Scientific Facility Program; CAS un- 
der Contracts Nos. KJCX2-YW-N29, KJCX2-YW-N45; 
100 Talents Program of CAS; Istituto Nazionale di Fisica 
Nucleare, Italy; Russian Foundation for Basic Research 
under Contracts Nos. 08-02-92221, 08-02-92200-NSFC- 
a; Siberian Branch of Russian Academy of Science, joint 
project No 32 with CAS; U. S. Department of Energy 
under Contracts Nos. DE-FG02-04ER41291, DE-FG02- 
91ER40682, DE-FG02-94ER40823; University of Gronin- 
gen (RuG) and the Hclmholtzzentrum fuer Schwerionen- 
forschung GmbH (GSI), Darmstadt; WCU Program of 
National Research Foundation of Korea under Contract 
No. R32-2008-000-10155-0 



[1] G.T. Bodwin, E. Braaten and G.P. Lepage, Phys. Rev. 

D 51, 1125 (1995); H.-W. Huang and K.-T. Chao, Phys. 

Rev. D 54, 6850 (1996); A. Petrelli, Phys. Lett. B 380, 

159 (1996); J. Bolz, P. Kroll and G. A. Schuler, Eur. 

Phys. J. C 2, 705 (1998); S.H.M. Wong , Eur. Phys. J. C 

14, 643 (2000). 
[2] C. Amsler and F.E. Close, Phys. Rev. D 53, 295 (1996). 
[3] M. Ablikim et al. (BESIII Collaboration), Nucl. Instrum. 

Meth. A 614, 345 (2010). 
[4] J.Z. Bai et al. (BES Collaboration), Nucl. Instrum. Meth. 

A 344, 319 (1994); J.Z. Bai et al. (BES Collaboration), 

Nucl. Instrum. Meth. A 627, 319 (2001). 
[5] K. Nakamura et al. (Particle Data Group), J. Phys. G 



37, 075021 (2010). 

[6] S. Agostinelli et al. (geant4 Collaboration), Nucl. In- 
strum. Meth. A 506, 250 (2003). 

[7] S. Jadach, B.F.L. Ward and Z. Was, Comput. Phys. 
Commun. 130, 260 (2000); S. Jadach, B.F.L. Ward and 
Z. Was Phys. Rev. D 63, 113009 (2001). 

[8] D.J. Lange et al, Nucl. Instrum. Meth. A 462, 1 (2001). 

[9] M. Ablikim et al. (BESIII Collaboration), Phys. Rev. 
D 81, 052005 (2010). 
[10] N. Berger, K. Zhu et. al., Trigger efficiencies at BESIII, 
Chinese Physics C 34, 1779 (2010) 



